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Abstract 

The coupHng between Ricci scalar curvature and the baryon num- 
ber current dynamicaUy breaks CPT in an expanding universe and 
leads to baryon asymmetry. We study the effect of time dependence of 
equation of state parameter of the FRW universe on this asymmetry. 



1 Introduction 

The origin of the difference between the number density of baryons and 
anti-baryons is still an open problem in particle physics and cosmology. The 
measurements of cosmic microwave background [1], the absence of 7 ray 
emission from matter- antimatter annihilation [2| and the theory of Big Bang 
nucleosynthesis [3j indicate that there is more matter than antimatter in the 
universe. Observational results yield that the ratio of the baryon number to 
entropy density is approximately n^/s ~ 10^^'^. In [4], it was pointed out 
that a baryon-generating interaction must satisfy three conditions in order 
to produce baryons and antibaryons at different rates: (1) baryon number 
violation; (2) C and CP violation; (3) departure from thermal equilibrium. 

In [5], by introducing a scalar field coupled to baryon number current 
it was suggested that the baryon asymmetry may be generated in ther- 
mal equilibrium while the CPT invariance is dynamically violated. Simi- 
larly , in [6], by introducing an interaction between Ricci scalar curvature 
and baryon number current which dynamically violates CPT symmetry in 
expanding Friedman Robertson Walker (FRW) universe, a mechanism for 
baryon asymmetry was proposed. The proposed interaction shifts the en- 
ergy of a baryon relative to that of an antibaryon, giving rise to a non-zero 
baryon number density in thermal equilibrium. The model suggested in 
O |6] was the subject of several studies on gravitational baryogenesis and 
leptogenesis in different models of cosmology 
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But, in [6], the problem was restricted to the cases that the equation 
of state parameter of the universe, u, is a constant, and the role of time 
dependence of to in baryogenesis was neglected. As a consequence, in this 
scenario, the baryon number asymmetry cannot be directly generated in 
radiation dominated epoch. But in [8j, in the framework of modified the- 
ories of gravity, following the method of [6], it was shown that the baryon 
asymmetry may be generated even in the radiation dominated era. 

In this paper, like [6], we assume that the universe is filled with perfect 
fluids such as the scalar infiaton field and radiation. The time dependence 
of uj is due to the fact that: (1) these components have different equation of 
state parameters, (2) they interact with each other and, (3) they may have 
time dependent equation of state parameters. We will study the effect of 
each of these subjects on time derivative of the Ricci scalar and consequently 
on baryogenesis. We will elucidate our discussion through some examples. 

Natural units h = c = ks = I are used throughout the paper . 



2 The role of u in gravitational baryogenesis 

The mechanism proposed in |6J to generate baryon asymmetry is based on 
the introduction of a CP violating interaction between the derivatives of the 
Ricci scalar and the baryon number current, J^, given by 

/ d'x^id,R)J>^, (1) 



M2 

where is a cutoff characterizing the scale of the energy in the effective 
theory and e = ±1. To obtain the chemical potential for baryon (hb) and 
antibaryons (/x^) for spatially constant R, we use 

^id,R)J^ = ^R{nB-ns), (2) 

where ub and are the baryon and antibaryon number densities respec- 
tively. Therefore the energy of a baryon is shifted by an amount of 
relative to an antibaryon. Thereby the interaction ([T]) dynamically violates 
CPT. We can assign a chemical potential to baryons: ^b = —fJ-B = ~^'m^- 
So, in thermal equilibrium there will be a nonzero baryon number density 
given by : 

9bT^ f 2^J'B , ^^^B.z 



nb = nB-n^ = -^[TT — + { — ) y (3) 

where gi, ~ 0{1) is the number of internal degrees of freedom of baryons. 
The entropy density of the universe is given by s = ^g-g'^T^, where Qs — 106 
indicates the total degrees of freedom for relativistic particles contributing to 
the entropy of the universe ^ . In the expanding universe the baryon num- 
ber violation decouples at a temperature denoted by To and a net baryon 
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asymmetry remains. The ratio ^ in the hmit T ^ mi, {nib indicates the 
baryon mass), and T S> ^j, is then: 

nb _ ^59b R I /.X 
s ~ ^47r25,M2T'^^" ^ ' 

Note that in different models we may have i? < as weU as R > 0, therefore 
introduction of e gives us the possibihty to choose the appropriate sign for 
rib- 

The geometry of the universe is described by the spatiahy flat FRW 
metric 

ds^ = dt^ -a^{t){dx^ + dy^ + dz^), (5) 
where a{t) is the scale factor. The Hubble parameter, H = ^, satisfies 

3 ^ 

H = -4ttG{P + p). (6) 

P and p are the pressure and energy density. We assume that the universe, 
filled with perfect fiuids, satisfies the effective equation of state: P = top. 
The equation of state parameter, can be expressed in terms of the Hubble 
parameter: uj = —1 — j^. The Ricci scalar curvature is given by 

= -3^2(1 -3w). (7) 

From Eq. ([6]), it follows that 

^ = ^(1 + a;)(l - 3^)pi + -^pdj. (8) 

Mp ~ 1.22 X 10^^ Gev is the Planck mass. If a; = 0, Eq. 1^ reduces to the 
result obtained in [6]. 

In the following we continue our study with a universe dominated by two 
perfect fluids with equation of states P^ = 'jdPd and Pm = Impm interacting 
with each other, through the source term Tipd + T2Pm'- 



Pd + 3H{pd + Pd) = Eipd + T2Pn 
Pm + 3i?(/9m + Pm) = -^iPd " ^2pm (9) 

Although these components don't satisfy the conservation equation solely, 
but 

p + m{i+uj)p = {). (10) 
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Note that p ~ pm + Pd and P ~ Pm + Pd- can be written in terms of the 
ratio of energy densities, r = — , 

w = — ; • (11) 

1 + r ^ ' 

From 

Pm Pd 

r = r — , (12) 

Pd Pd 

and Eq. ([9]) one can determine the behavior of r with respect to the comoving 
time 

r = -{r + l){Ti + Tsr) + 3ifr(7d - 7^), (13) 
or by suppressing 7^ 

r = -(r + l)(ri + Tsr) + ?,Hr{l + r){uj - 7„). (14) 
p and r may be related through 

p r + (r + l)(ri + T2r) , , . . , . 

Combining Eqs. (|lip and (|12p results in 

. _ id + rjrn _ (7m - 7rf)(ri + ^2^) _ 3 (7m - ld)'^Hr 
^~ l + r 1 + r (l + r)2 ■ ^ ' 

The first and second terms show the effects of time dependence of 7's, and in- 
teraction of fluid components on Co, respectively. Note that even for constant 
7's and in the absence of interactions, as the third term of (jl6p indicates, 
u) varies with time. This is due to the fact that the universe is assumed 
to be composed of components with different equation of state parameters. 
Putting Eq. (fT6l) into Eq. ([8]) gives 



A _ \/3 / (1 - 37m) (1 + 7m)r^ + 2(1 - 7d - 7m - '?>ldlm)r 

(1 + r)2 

(7d + l)(-37rf + l)\ 3 3 ( id + rim \ 



[1 +r)2 j"^ \ 1 + r 



3 (El + rr2)(7m - 7d) 3\/3 (7m - 7^)^?- | 



M2 l + r M3 (1 + r) 

This equation can be rewritten in terms of uj 

n ^n_L Q ^ I _L (7m - uj)id + {uJ- Jdh'm 
''=W^^'+^^^'- '^^P' + Ml ^i-^, P - 

3 3\/^ 3 

—2 ((7m - + (W - 7d)r2) p - -J7^{^ - lm){ld " . (18) 
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If we neglect to, as mentioned before, only the first term remains: R = 
-^(1 + uj){l — 3uj)p2 which is zero at w = | and at a; = —1. But by taking 

UJ into account we may have baryon asymmetry at w = | and lo = —1. 
Although the asymmetry generated during inflation will be diluted away. 

It is worth to study what happens when one of the fluid components 
corresponds to radiation (e.g. produced after the inflation epoch). To do 
so, we take Am = ■j- this case Eq. ([T7|) reduces to 

• ^ V3 (l-37rf)(l+7rf) I 1 (l-37rf)(ri + r2r) 3 id_ 

M3 (1 + r) M2 1 + r ^+M2l + r^' 

(19) 

In general F's may also be time dependent [10], e.g., one can consider Fi = 
XiH and r2 = X2H , where Ai, A2 € [TT]. Depending on the model under 
consideration, the third term of (jl9p . including the time derivative of 
may be simplified in terms of other parameters of the model. For example 
consider a massive scalar field of mass m, with a time dependent equation 
of state parameter interacting with radiation. The energy density, pd, and 
the pressure, Pd, of the scalar field satisfy 

■'■ j,2 I l 2/2 

Pd = (p 

Pd = l^'-l^n'^'- (20) 
The time dependent equation of state parameter of the scalar field is then 

_ i</>2 - Im^^' ^^^^ 



The scalar field interacts with another component (radiation) characterized 
by the equation of state parameter 7^ = |: 

Pd + 3i?(/0d + Pd) = Tipd + T2Pn 

pn + ^Hpn = -^lPd-^2P7^ (22) 

The subscript TZ denotes the component with 7^ = . To determine ''jd in 
Eq. (fT9|l . let us define n = (1 — 7d)^- Then using 



ii = mpd^l-ll (23) 
which was verified in [12], we are led to 



id = -2m^l - 7i + (1 - 7d)(ri + rr2) - 3F(1 - Yd)- (24) 
By supposing 

ri = 3ArfF, ^2 = 3A7^F, Ad,A7^G3^^, (25) 
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we arrive at 



id = -2m^l - 72 + 3/7(1 - 7d)(Ad + Xnr -Id -I)- (26) 
Finally we deduce 



. 6mVl-Td , 2^(-l-7rf + Ad + rA7^) 3 



If the potential is negligible with respect to the kinetic energy, (j) ^ m 
we have 7^ ~ 1 and 7^ ~ 0. In this case R=^ (^^±^+7^) Pd- 
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When the scalar field dominates, i.e., in the limit r — > 0, = — ffs- y 1 — l'dPd+ 

^(-l-7d + A,)4. 

In the limit r ^ 00, i.e, when dominates 

^=1^Vp|- (28) 

p 

So, even in the radiation dominated era, we can have i? 7^ if X-ji 7^ 0. 
For radiation component, the energy density is related to the equilibrium 
temperature via p| 

pn = enT\ (29) 

2 

where eTi = ^Qi, and counts the total number of effectively massless de- 
grees of freedom. Its magnitude is g^, Qs- Therefore the baryon asymmetry 
in terms of temperature can be determined from (j28p . and ()29p : 

T-2.5»a.-j||j. (30) 

By defining = aMp we find that ^ is of order 

5 

(31) 

When 7's are time dependent, in general, derivation of the exact solutions 
of Eq. ()22p is not possible. Thus explaining the precise behavior of the fluid 
components in terms of temperature, except in special situations like (jSOp . 
is not straightforward. But when 7's are constant, it may be possible to 
determine baryogenesis in terms of decoupling temperature. We will study 
this situation through two examples. 

As a first example consider 7^ = 77^ = 1/3. We assume that 7^ is 
a constant satisfying 7^ > | corresponding to a non-thermal component 
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decreasing more rapidly than radiation p]. If there is no energy exchange, 
i.e, Ti = = 0: 

Pd + 3{jd + l)Hpd = 

PTi + AHpTi = 0. (32) 

In this case 

^ = i|(l^Mfi±l^),i (33) 

(l + r) ^ ' 

This result may reproduced by taking 7,^ = 1 (which results in 7^ = 0) and 
-^Te = Arf = in (p7|) . r satisfies 

r = Hr{3^d - 1), (34) 

hence, r is an increasing function of comoving time, i.e., pd decreases faster 
than the radiation component. R is also an increasing function of time. To 
see this let us consider 

^ = ^(1 + 7d)(37d - 1) (9(7d + 1) + 2r(37d + 5)) ^^^^ (35) 

which was obtained using Eqs. (fT5|) and ([33]) . Hence "Yd > h leads to i? > 0. 



Eq. ([5^ implies pn (x a ^, hence the temperature redshifts as T 



oc a 



-1 



this relation and Eq. ()32p imply /j^ oc T^(^+'>'<*). If at a temperature denoted 
by r = Trd, pn = Pd, we can write 

Pd = (--rTrd ( j , (36) 



which results in 



and 



(38) 



by putting Eqs. ^ and in Eq. ([M]) and subsequently in Eq. (g]) at 

T = To, the baryon asymmetry may be determined: 

In the limit r — > (related to Td ^ Trd), 

fj 9(l+7d) 
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which is the same as the result obtained in where the effect of time 
dependence of uj was ignored, indeed in this hmit lo does not contribute in 
baryon asymmetry. In the hmit r ^ 1 which corresponds to pd <^ pn we 
have 

E. g. if we define Td = /3Tr£), /? ^ 1, and as before = aMp, then the 
order of magnitude of ^ is obtained as 

As a second example consider a model, with constant 7^, and 7^ = g- 
Also assume that at least one of the constant F's are not zero. Hence 

^ (l-37d)(l+7d) J 1 (l-37d)(ri + r2r) _^ 

^ = OTO ' TTT ^' (^^^ 

may be used to determine In this situation although Eq. ([^U]) holds but 
Eq. (f36]) is no longer valid. The temperature does not behave as T oc 
and, instead, we have T = f{a) [S], [13]. In some special cases, R can be 
easily determined in terms of decoupling temperature. For example, at r ~ 1 
(when UJ ~ ^-^^-^) we have 

P ^ 2enT^ (44) 

and 

T - f (^P^- + ^^'^S - 5|(3.. - i)(r, + r.)Tg) . (45, 

Note that in the above example Tq = T^o. 

Also For r — > and r — > 00, corresponding to a; ~ 7^ and a; ~ 7^ 
respectively, Eq. (fTB]) leads to 

[ (7d - 7m)r2, r ^ 00 

Therefore for 7^ = |, 

. _ / ^(1 + 7.)(1 - 37.)4 - ^i^P., - - 

E.g., if one suppose that the reheating process after the infiation is due to 
the decay of a scalar field (the infiaton field: 0) to (ultra)relativistic species, 
via coherent oscillation [9], then 

PR + AHpn = T^p^, (48) 
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where pn is the energy density in the relativistic decay products. If 
is the vacuum energy of the scalar field at the beginning of the oscillation, 
from t ~ until t ~ T^^, (f) particles, which behave like non-relativistic 
particles (i.e., 7,/, = 0), dominate the energy density [9j. During this time, 
a{t) oc 1 3 , p (X and 

Although pti oc and Eq. ()29p is still valid, but due to the interaction 
term, pTi oc does not hold. Then T (x is not valid. In this case one 
can show that the temperature can be approximated by [13] 



TRn{^f. (50) 



Hence at the time of decoupling 



~ i^^^^ + ivf ^^^^^ ) ■ ^^^^ 

To derive the above equation we have multiplied our result by the dilution 
factor [6j. The first term is the same as the result obtained in [6], 

where io was ignored. If we take Ta. ~ , and as before define A^~^ — cxAlp 
and Td = PTru, then we are able to compare the order of magnitudes of 



and 



terms contributing in (I52p : 




' '''^TnoTh]~i(^]\ (54) 



When /3 ^ 1 the second term in (I52p which comprises the effect of decay 
width is negligible. 

At the end let us note that if to were time independent, ^ would be 
vanish in radiation dominated era. But by taking uj into account, it is easy 
to see that if r2 / 0, in the limit p-ji ^ p^ we have R = ^^-^^"^^ pn^ 
therefore 

nb 9b (37d - l)r2 



Thereby ^ ~ ^''S'^^lJf - 
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